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基于 济 固 耦合 的 条 汽 发 生 器 传 热 管 结构 应 力 分 析 


SUE BE 史 强 * Kbm 唐 晖 XIIA 


《生态 环境 部 核 与 辐射 安全 中 心 ”北京 100082) 


要 ”蒸汽 发 生 器 传 热管 是 压 水 堆 核 电站 一 回路 的 关键 压力 边界 , 正常 运行 及 事故 条 件 下 的 系统 运行 参数 ， 
接 影 响 着 传 热管 结构 的 完整 性 。 为 进一步 对 传 热管 的 结构 响应 特性 进行 研究 ， 建 立 了 两 种 主流 排列 方式 的 
管 局 部 模 
H ] 


传 热 型 ， 并 分 别 进行 了 正常 运行 和 典型 事故 参数 下 的 流 固 耦合 计算 和 对 比分 析 ， 获 得 了 一 、 二 次 侧 
流体 共同 作用 下 的 传 热管 关键 部 位 的 载荷 变化 规律 。 结 果 表明 ， 在 相同 的 流动 条 件 下 ， 传 热管 又 排 模型 的 应 
力 和 变形 均 大 于 顺 排 模型 传 热管 根部 横 截面 的 等 效应 力 沿 壁 厚 方向 先 减 小 再 增 大 ， 中 部 模 截 面 的 等 效应 力 

沿 壁 厚 方向 呈现 近似 线性 降低 ， 在 传 热管 根部 区 域 ， 内 、 外 壁面 的 等 效应 力 沿 轴 向 高 度 先 迅 速 降低 而 后 快速 
升 高 ， 并 形成 应 力 平台 


E 
[zi 


nk 
C 


关键 词 KARE: 传 热管 ， 流 固 耦 合 ， 应 力 分 析 
中 图 分 类 号 TL364 
DOI: 


Stress Analysis of Heat Transfer Tube Structure in Steam Generator Based on Fluid 


Structure Coupling Method 


Guo Chao Wei Chao Shi Qiang* Zhu Guangyu Tang Hui Liu Qiaofeng 
(Nuclear and Radiation Safety Center, Beijing, 100082, China) 


Abstract [Background] The heat transfer tube of steam generator is critical pressure boundary in the primary 
circuit of the pressurized water reactor nuclear power plant. The system operating parameters during normal operation 
and accidents directly affect the integrity of the heat transfer tube structure. [Purpose] This study aims to obtain the 
structural response characteristics of heat transfer tube models with different arrangements under normal operation 
and accident system operating parameters. [Methods] Local models of heat transfer tube with two mainstream 
arrangements were established, and fluid structure coupling calculations and comparative analysis were carried out 
under normal operation and typical accident parameters, respectively. [Results] The maximum stress and deformation 
values of heat transfer tubes under normal operation and extreme accident parameters were obtained. The load 
distribution patterns of typical axial and radial positions of heat transfer tubes were obtained. [Conclusions] The 
results show that under the same flow conditions, the stress and deformation of the heat transfer tube cross 
arrangement model are greater than those of the straight arrangement model. The equivalent stress at the root 
cross-section of the heat transfer tube decreases first and then increases along the wall thickness direction, while the 
equivalent stress at the middle cross-section shows an approximate linear decrease along the wall thickness direction. 
In the root area of the heat transfer tube, the equivalent stress on the inner and outer walls first rapidly decreases and 
then rapidly increases along the axial height, forming a stress plateau. 
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过 程 中 由 于 
蒸汽 发 生 器 是 压 水 推 


要 影响 。 
要 位 置 [1 


增 大 的 事故 ， 尤 其 是 压 差 达到 


发 更 加 严重 的 放射 性 泄露 后 果 辐 。 因 此 ， 


F 温差、 压 差 和 各 类 化 学 腐蚀 等 作 月 


包 广 的 关键 设备 之 一 ， 是 阻止 放射 性 
上 月， 不 可 避免 的 出 现 老 化 降 质 现象 ， 会 对 ! 
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事故 的 四 


2, ZjidET XU PER 


回路 压力 边界 最 薄弱 的 环节 
核电 站 燕 汽 发 生 器 正常 运 
JIZ 15 MPa 和 6MPa， 较 大 的 压 差 可 能 会 引发 传 热管 变形 ， 现 有 研究 表明 导致 蒸汽 发 生 器 一 、 
11MPa 以 上 时 ， 极 易 引 发 蒸汽 发 生 器 传 热 


分 之 一 。 压 水 堆 


研究 一 、 二 次 侧 流 


汽 发 生 器 的 优化 设计 具有 十 分 重要 的 意义 。 


上 述 过 程 主要 涉及 一 次 侧 流体 与 传 热管 内 壁 、 
些 学 者 利用 流 
通过 对 蒸汽 发 生 器 管 进 行 双 向 流 固 耦 合 
体 与 换 热 管 
布 ， 郑 陆 松 等 四 
EE 


最 高 应 力 值 及 某 些 路 径 的 应 力 变化 规律 ，Sun Baozhi 等 


玉 爽 等 


传 热管 儿 


回 耦 合 方法 开展 了 传 热管 


及 文 撑 板 的 耦合 
进行 了 蒸汽 发 生 器 一 、 
H vi 2 EHS HJ 77 4G] RE 


静态 结构 分 析 ， 计 算 了 热 应 力 和 力 载荷 下 的 机 械 应 力 以 及 两 种 载荷 
常 运行 参数 下 的 传 热管 位 移 和 应 力 分 布 特点 。 

现 有 研究 主要 针对 特定 排列 方式 的 传 热管 〈 顺 排 为 主 ) j 
研究 未 同时 考虑 一 、 二 次 侧 流 体 的 共同 作用 。 本 文 基于 某 三 
和 又 排 两 种 主要 排列 方式 的 传 热管 模型 ，j 
算 和 对 比分 析 ， 获 得 了 一 、 二 次 侧 流体 共 


计 及 运行 提供 参考 。 
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蒸汽 发 生 器 传 热管 物理 模型 


EKHE 


合流 动 换 热 进行 三 维 
二 回路 耦合 数值 模拟 ， 对 
换 热 器 换 热 管 及 壳 体 相连 接 的 管 板 ; 
等 090 通 过 流 


的 力学 响应 研究 了 


[人 


， 获 得 了 传 热管 的 运动 规律 ; 
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物质 泄漏 的 重要 屏障 ， 蒸 汽 发 生 器 在 运行 
安全 运行 产生 重 
和 故 在 核 动力 厂 设备 事故 中 居 
二 次 侧 流体 平均 压 
二 次 侧 压 差 
管 破裂 (SGTR) B4， 进 而 可 能 引 
管 结构 响应 特性 ， 对 蒸 
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， 其 传 热管 断裂 
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场 及 压 差 作用 下 的 传 热 


壁 与 二 次 侧 流体 之 间 的 双向 流 固 耦 合 问题 
Æ, W Kuehlert 等 四 在 考虑 二 次 侧 流体 的 作用 下 
孙 宝 芝 等 中 利用 CFX 对 一 、 二 次 侧 流 


TUBE. f 


究 了 换 热 管 在 流体 压力 场 载荷 作用 下 的 应 力 分 
车 热 管 进行 稳 态 热 分 析 和 热 应 力 分 析 ; 陈 
行 了 应 力 分 析 ， 给 出 了 管 板 
行 了 传 热管 和 支撑 板 的 热 分 析 和 
合 应 力 等 ， 赵 颖 杰 等 0 研究 了 正 
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核电 站 燕 汽 发 生 器 的 类 型 主 


列 管 式 直 流 蒸 汽 发 生 器 等 ,其 中 立 式 U 型 管 
环 蒸汽 发 生 器 一 般 | 
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型 管束 、 


3n 


TE 


龙 一 号 机 型 


管 板 、 防 震 条 以 及 管束 支 
角形 ) 两 种 方式 (图 1), 如 M310 机 型 
4 采用 又 排 排列 ，CPR1000 机 型 


HALA U 型 管 自然 循环 蒸汽 发 生 器 、 著 式 
目 然 循环 薰 汽 发 生 器 在 国内 应 有 


上 下 封 头 、 


p. 
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4 的 ZH-65 型 采用 又 排 排列 等 。 


管束 区 域 、 
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55/19B # 


Jf] 55/19B 型 采用 顺 排 排列 ，AP1000 机 型 
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a) 顺 排 排列 


“先进 核能 技术 设计 


与 安全 教育 部 


并 分 别 进行 了 正常 运行 和 典型 事故 参数 下 的 流 回 耦合 
同 作 用 下 的 传 热管 关键 部 位 载荷 变化 规 得 


汽水 分 离 装置 和 简体 等 结构 组 成 ， 其 中 管束 区 域 具体 
热管 的 排列 方式 主要 有 顺 排 〈 正 方形 ) 和 又 排 〈 三 
型 燕 汽 发 生 器 传 热 管 采用 顺 排 排列 , EPR 机 型 


行 正常 运行 参数 下 的 力学 响应 分 析 ， 且 部 分 
动 核 电站 蒸汽 发 生 器 原型 ， 建 立 了 顺 排 
合 数值 模拟 计 
EE， 可 为 蒸汽 发 生 器 的 设 


代 非 能 


自然 循环 蒸汽 发 生 器 和 
i£. AU 型 管 自然 循 
包括 有 U 


最 为 
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b) 又 排 排列 


图 1 传 热管 排列 方式 
Fig. 1 Arrangement of heat transfer tubes 
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郭 超 等 ， 基于 流 固 耦合 的 蒸汽 发 生 器 传 热管 结构 应 力 分 析 


由 于 传 热管 结构 复杂 、 管 束 细 长 且 数 量 可 达 上 万 根 ， 对 传 热管 束 的 分 析 的 一 般 选 取 单 元 管 和 周边 相 邻 
管束 轮廓 ， 以 及 相应 的 一 、 二 次 侧 流体 域 来 进行 物理 模型 的 简化 ， 如 图 2 虚线 框 内 所 示 。 考 虑 传 热管 的 根部 
位 置 是 蒸汽 发 生 器 事故 的 多 发 区 023， 本 研究 选取 了 蒸汽 发 生 器 热 端的 管 板 连接 处 ， 沿 传 热管 高 度 Im《〈 即 相 
邻 两 块 管束 支撑 板 间 距 ) 的 局 部 结构 为 主要 研究 对 象 ， 模 型 的 几何 参数 参考 某 三 代 非 能 动 核电 站 蒸汽 发 生 
器 原型 ， 蒸 汽 发 生 器 的 U 型 传 热管 内 径 为 15.44mm， 壁 厚 为 1.02 mm， 管 间距 为 24.89mm 。 传 热管 材料 为 
Inconel 690 合 金 ， 其 弹性 模 量 为 211GPa， 届 服 强度 为 330MPa， 泊 松 比 为 0.289， 密 度 为 8190kg/m3。 建 立 的 
顺 排 和 又 排 单元 管 计 算 模型 如 图 2 所 示 ， 两 种 模型 中 的 一 次 侧 和 传 热管 的 几何 尺寸 相同 ， 二 次 侧 的 流动 截面 
不 同 ， 顺 排 模型 截面 积 为 1318.12mm2， 叉 排 模型 为 889.60mm2。 


Outer wall profile of 
adjacent tubes 


Secondary 
side fluid 


Heat transfer 


Primary 
side fluid 
图 2 传 热管 顺 排 〈 左 ) 和 又 排 〈 右 ) 单元 管 模型 


Fig.2 Unit tube straight arrangement model (left) and cross arrangement model (right) 


2 ”模型 网 格 验 证 及 计算 工 况 
2.1 ”边界 条 件 


对 于 单元 管 模型 的 固体 域 ， 考 虑 管束 下 端 一 般 固 定 在 直径 Sm, EE REOR BI SS LE. H 
沿 高 度 等 间距 I 1m) 有 管束 支撑 板 进 行 约 束 ， 因 此 对 所 研究 的 管束 段 模型 的 上 、 下 端面 直接 设置 位 移 
固定 约束 。 将 传 热管 外 壁 和 二 次 侧 的 接触 面 ， 以 及 内 壁 和 一 次 侧 的 接触 面 设置 为 流 固 耦合 作用 交界 面 ， 以 
实现 一 、 二 次 侧 对 管束 的 共同 作用 。 将 模型 外 边界 的 相 邻 传 热管 轮廓 面 设 置 为 壁面 边界 ， 其 余 的 模型 外 边 
界 均 设 置 为 对 称 边界 。 对 于 模型 的 流体 域 ， 在 计算 过 程 中 对 一 、 二 次 侧 进口 采用 流速 边界 条 件 ， 一 、 二 次 
侧 出 口 采用 压力 边界 条 件 ， 参 数值 依据 具体 的 计算 工 况 设置 。 


2.2 ”网 格 无 关 性 验证 


建立 了 02m 的 顺 排 传 热 单元 管 模 型 进行 了 网 格 无 关 性 验证 。 通 过 对 模型 进行 合理 切 分 后 ， 再 进行 网 
格 划分 以 保证 网 格 质量 。 对 一 、 二 次 侧 流体 域 分 别 进行 网 格 平均 尺寸 为 Imm, 1.5mm, 2mm. 3mm, 4mm 
的 网 格 划分 ， 其 中 ， 对 于 流体 域 为 Imnm 和 1.5mm 网 格 尺寸 的 模型 ， 其 相应 的 传 热管 在 厚度 方向 设置 三 层 
网 格 ; 其 余 网 格 尺寸 的 模型 对 应 传 热 管 厚度 方向 尺寸 为 两 层 网 格 。 对 上 述 5 组 网 格 尺寸 模型 分 别 进行 了 蒸 
汽 发 生 器 正常 运行 参数 的 流 固 耦合 计算 ， 计 算 完 成 后 提取 传 热管 根部 截面 的 壁面 厚度 中 点 沿 轴 向 的 Von 
Mises 等 效应 力 ， 结 果 如 图 3 所 示 。 当 流体 区 域 网 格 尺 寸 达到 1.Smm， 传 热管 厚度 方向 为 三 层 网 格 时 ， 进 
一 步 细 化 网 格 对 计算 结果 几乎 无 影响 ， 结 果 已 趋 于 稳定 ， 兼 顾 计算 精度 与 效率 ， 最 终 选 定 模型 平均 网 格 尺 
十 为 1.Smm。 最 终 所 建立 的 蒸汽 发 生 器 传 热管 根部 单元 管 模 型 及 网 格 划 分 如 图 4 和 图 5 所 示 。 
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Flow direction (axial direction) / mm 


图 3 不 同 网 格 


尺寸 模型 沿 轴 向 高 度 的 等 效应 力 


Fig.3 Von-Mises stress of models with different grid sizes along the axial height 
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a) 流体 网 格 


图 4 顺 排 单元 管 模型 及 网 格 划 分 
Fig.4 Model and Grid Division of Straight Unit Tube 


计算 工 况 


b) 固体 网 格 


输入 参数 ， 该 时 刻 主 泵 已 


A — 


FIZ, 
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c) 流体 网 格 


图 5 又 排 单元 管 模型 及 网 格 蕊 
Fig.5 Model and Grid Division of Cross Unit Tube 


d) 固体 网 格 
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[ 况 ， 选 取 了 电厂 满 功 率 运 行 的 一 、 二 次 侧 流 速 和 压力 参数 作为 计算 输入 参数 。 对 二 
故 工 况 ， 选 取 了 可 导致 一 、 二 次 侧 压 差 显 著 增 大 的 两 类 典型 始 发 事故 ， 即 属于 一 次 侧 热 输出 增加 类 的 主 
汽 管 道 断裂 事故 (Main Steam Line Break, MSLB )， 以 及 属于 二 次 侧 系 统 热 移 出 能 力 减 小 类 的 
(Loss of Feed Water LOFW )。 对 于 MSLB 工 况 ， 根 据 该 电 
侧 压 差 最 大 时 刻 的 系统 参数 作为 计算 输入 参数 ， 该 时 刻 主 泵 
尚未 停 堆 ， 一 次 侧 压力 仍 处 于 高 位 ; 而 该 时 刻 3 


最 终 安全 分 析 报 告 ， 选 取 破 损 
尚未 关闭 ， 一 次 侧 仍 处 于 正常 运行 流速 ， 由 
蒸汽 隔离 阀 已 关闭 ， 二 次 侧 流速 处 于 较 低 水 平 ， 
力也 由 于 破 口 喷 放 而 降低 。 对 于 LOFW 事故 ,选取 含 PRHR 环 路 的 一 、 二 次 侧 压 差 最 大 时 刻 的 系统 参数 人 
管束 一 次 侧 进口 流速 显著 降低 ， 由 于 该 时 刻 处 于 PRHR 投入 后 传 热 
在 逐渐 降低 的 阶段 ， 一 次 侧 压 力 则 处 于 高 位 ， 而 由 于 该 时 刻 汽轮机 已 停机 ， 二 次 侧 流 速 则 基本 为 零 ， 
则 压力 也 处 于 整个 事故 阶段 的 较 低 水 平 。 正 常 运行 和 事故 ] 


[ 况 的 具体 计算 参数 如 下 : 
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MSLB condition 
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表 1 正常 运行 和 事故 工 况 的 计算 参数 
Table 1 Calculation parameters for normal operation and accident conditions 
ds 一 次 侧 流速 一 次 侧 压力 二 次 侧 流 速 二 次 侧 压 力 

工 况 : ; : i : 

condition Primary flow Primary side Secondary side Secondary side 
velocity pressure flow velocity pressure 
Ti- 
Se 0.5 m/s 15.50 MPa 0.1688 m/s 5.60 MPa 
Normal operation 
MSEE TUR 0.5 m/s 15.20 MPa 0.0314 m/s 5.03 MPa 


二 次 侧 压 


丧失 给 水 事故 
路 一 、 二 次 


SpA: dE TARS RU ART UR IE SR TRAVES AIL] 2T 
LOFW 工 况 近似 为 0〈 取 
LOFW condition 0:0213 ms AMP amimia OP 


3 ”计算 结果 分 析 
3.1 ”整体 结构 响应 


对 正常 运 
变形 量 最 大 ; 
MSLB 和 LOFW 工 况 参数 下 的 最 大 应 力 和 变形 值 对 比 。 


表 2 正常 和 事故 工 况 传 热管 计算 结果 


Table 2 Calculation results of heat transfer tubes under normal and accident conditions 


Z11 fl MSLB. LOFW 参数 下 的 传 热管 流 回 耦合 计算 结果 表明 ， 传 热管 结构 的 根部 区 域 应 力 及 
结构 的 整体 响应 以 传 热管 中 部 横 截 面 基 本 呈 对 称 分 布 。 表 2 为 顺 排 和 又 排 结构 在 正常 运 
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大 等 效应 力 最 大 变形 量 
Maximum Von-Mises stress / MPa Maximum deformation / mm 
顺 排 又 排 偏差 顺 排 X HE 偏差 
Straight Cross Deviation Straight Cross Deviation 
em 110 112 1.8% 0.28x102 029x102 3.6% 
Normal operation 
MSLB 工 况 3 3 E " 
MSLB condition 113 115 1.896 0.29x10 0.31x10 6.9% 
LOFW Li 3 z 3 7 
LOFW condition 143 146 2.1% 0.37x10 0.40x10 8.1% 
从 表 可 以 看 出 , 在 正常 运行 和 事故 工 况 下 , 顺 排 模 型 的 最 大 应 力 及 变形 量 均 小 幅 大 于 顺 排 模 型 。 LOFW 
事故 工 况 下 的 结构 响应 更 加 显著 ， 其 顺 排 模型 最 大 应 力 和 变形 比 正常 运行 分 别 高 出 30.0% 和 30.3%; 又 排 


模型 比 正 常 运行 分 别 高 出 32.1% 和 37.9%。 


管内 壁面 ， 右 侧 为 外 壁面 。 图 a)、b) 和 c) 为 顺 排 模型 的 正常 运行 、MSLB fll LOFW 工 况 计算 结果 ; 
DMO 为 又 排 模 型 相应 计算 结果 ， 从 图 中 可 以 看 出 ， 传 热管 内 壁面 的 根部 区 域 应 力 最 大 。 


E | j | | j | 


Contour 2 
1.460e*008 
1.339e+008 
1.218e+008 
a) b) c) d) (e 
6 传 热管 根部 纵 截 面 应 力 分 布 云图 
Fig.6 Stress distribution of the longitudinal section at the root 
of the heat transfer tube 


1.097e+008 
9.760e+007 
8.550e+007 
7.340e+007 
6.130e+007 
4.920e+007 
3.710e+007 


2.500e+007 
[Pa] 


3.3 ”局 部 应 力 分 析 


根据 表 2 和 图 6， 顺 排 模 型 和 了 又 排 模 型 在 相同 工 况 参数 下 的 应 力 分 布 具 有 相似 性 ， 以 下 以 又 排 模型 
算 结 果 为 例 , 进一步 分 析 传 热管 的 局 部 载荷 特性 。 图 7 为 又 排 模 型 传 热 管 根部 和 中 部 横 截面 应 力 分 布 云 图 ， 
各 壁 厚 方向 应 力 先 降低 而 后 


根部 区 域 ， 内 壁面 应 力 最 大 ， 且 
XXXXXX-5 


从 图 7 的 a)、b) 和 c) 中 可 以 看 出 ， 在 传 热管 


图 6 为 传 热管 根部 区 域 纵 截 面 应 力 分 布 云图 ， 每 个 分 图 的 下 部 为 传 热 管 与 管 板 的 连接 处 ， 左 侧 为 传 热 


图 d). 


型 计 


Eh 技 术 20XX, XX: XXXXXX 


增 大 ， 从 图 7 的 D. e RED. 则 可 以 看 出 ， 传 热管 中 部 应 力 分 布 与 根部 不 同 ， 内 壁面 应 力 最 大 ， 沿 壁 厚 广 


向 应 力 逐 渐 降低 ， 外 壁面 应 力 最 低 。 


c) 根部 -LOFW 


Von Mises Stress Von Mises Stress Von Mises Stress 
;ontour 1 Contour 1 Contour 1 

1.1006«008 113364008 1.433e+008 
1.016e+008 1.046e+008 1.323e+008 
9.308e+007 9.590e+007 1.214e+008 
8.460e+007 8.719e+007 1.104e+008 
7.612e+007 7.847e+007 9.945e+007 
6.764e+007 6.976e+007 8.849e+007 
5.916e+007 6.104e+007 7.754e+007 
5.0686+007 5.233e+007 6.659e+007 
4.220e+007 4.361e+007 5.563e+007 
3.372e+007 3.490e+007 4.468e+007 

i picem 2.619e+007 3.373e+007 

[Pa] al 
立 
a) 根部 -正常 运行 b) 根部 -MSLB 

Von Mises Stre: Von Mi Stress Von Mises Stress 

Contour 2 Contour Contour 2 
8.025e+007 8.245e*007 1.035e«008 
7.840e«007 8.056e«007 1.012e«008 
7.655e+007 7.867e+007 9.887e+007 
7.470e+007 7.678e+007 9.653e+007 
7.286e+007 7.489e4007 9.419e«007 
7.101e+007 7.300e«007 9.185e+007 
6.916e+007 7.111e*007 B.951e*007 
6.314007 6.922e+007 8.717e+007 
6.546e+007 6.733e+007 8.483e+007 
6.361e+007 6.544e+007 8.249e+007 
6.176e+007 6.355e+007 8.015e+007 

[Pa] [Pa] [Pa] 


d) 中 部 -正常 运 


图 7 又 排 模型 根部 和 中 部 横 截面 应 力 分 布 云 


T e) 中 部 -MSLB 


D 中 部 -LOFW 


Fig.7 Stress distribution at the root and middle cross-sections of cross model 
为 进一步 展现 载荷 分 布 特性 ， 图 8 和 图 9 分 别 给 出 了 又 排 模型 传 热管 根部 和 中 部 横 截面 ， 沿 壁 厚 方向 
上 的 应 力 分 布 曲线 〈 径 向 坐标 原点 为 传 热管 中 心 )， 图 中 可 以 更 显著 的 看 出 前 文 所 述 的 应 力 分 布 特性 ， 传 热 
的 管 根部 在 三 种 工 况 下 的 最 低 值 较 内 壁面 应 力 降低 约 79.5MPa. 85.6MPa 和 107.6MPa; 较 外 壁面 应 力 分 别 


降低 约 31.8MPa、34.9MPa 和 45.1MPa; 事故 工 况 比 正常 运行 下 的 应 力 降低 和 升 高 过 程 更 迅速 ， 且 变化 幅 


值 也 更 大 。 传 热管 的 中 部 在 三 种 工 况 下 应 力 沿 壁 厚 呈 现 近似 线性 降低 ， 外 壁面 较 内 壁面 降低 约 16.8MPa、 


17.3MPa 和 21.8MPa。 


160 


—s— Normal 


140 -= MSLB 


—4— LOFW 


X 
N 
e 


=ñ 
o o 
O o 


Von-Mises stress / MPa 
Von-Mises stress / MPa 


7.6 7.8 8.0 8.2 8.4 8.6 8.8 


Wall thickness direction (radial direction) / mm 


160 
440 —s— Normal 
—e— MSLB 
120 —4— LOFW 
100 
80 
60 
40 
20 
0 一 一 
7.6 7.8 8.0 8.2 8.4 8.6 8.8 
Wall thickness direction (radial direction) / mm 


8 又 排 模型 传 热管 根部 沿 壁 厚 方向 的 应 力 分 布 图 9 又 排 模型 传 热管 中 部 沿 壁 厚 方向 的 应 力 分 布 
Fig.9 Stress distribution along the wall thickness 
direction at the root of heat transfer tube in cross model direction in the middle of heat transfer tube in cross model 


Fig.8 Stress distribution along the wall thickness 


10 和 图 11 221073 SCHEBEELPETAET P BETRURU PE ThE ER 


向 高 度 〈 即 流体 流动 方向 ) 上 的 应 力 分 布 。 


从 图 中 10 可 以 看 出 ， 内 壁面 根部 应 力 最 大 ， 且 沿 着 内 壁面 轴 向 高 度 ， 应力 首 先 迅速 减 小 再 迅速 增 大 并 形成 


现 与 内 壁面 类 似 的 变化 规律 ， 但 存在 以 下 差异 :一 方面 3 


局 部 应 力 峰 值 ， 而 后 呈现 应 力 平台 现象 ,平台 应 力 显 著 低 于 根部 应 力 。 在 图 11 P, WAIE, ME 


其 整体 应 力 水 平 人 


氏 于 内 壁面 ， 并 且 最 大 应 力 不 在 根 


部 ， 而 在 反弹 后 形成 的 应 力 峰 值 位 置 处 ; 另外 所 形成 的 平台 应 力 略 高 于 根部 应 力 。 图 10 中 三 种 工 况 下 的 平 
台 应 力 较 内 壁面 最 大 应 力 ， 分 别 降低 29.8MPa. 30.8MPa 和 60.4MPa; 图 11 中 三 种 工 况 下 的 平台 应 力 较 外 
壁面 最 大 应 力 ， 分 别 降 低 3.3MPa、3.4MPa 和 4.2MPa， 事 故 工 况 下 内 、 外 壁面 的 平台 应 力 与 应 力 最 大 值 的 
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: 基于 流 固 耦合 的 蒸汽 发 生 器 传 热管 结构 应 力 分 析 


差 值 大 于 正常 运行 。 
160 
—s— Normal 
140 —。— MSLB 
S —4— LOFW 
E 120 
$ 100 
$ 80 
三 
5 60 
- 
40 
0 5 10 15 20 25 30 
Flow direction (axial direction) / mm 
图 10 又 排 模型 内 壁面 轴 向 应 力 分 布 特性 
Fig.10 Characteristics of axial stress distribution on 
the inner wall of cross model 
4 结语 
本 文 基于 流 固 耦 合 方法 ,考虑 一 、 二 次 侧 流 体 的 


共同 作用 , 对 蒸汽 发 生 器 不 同 的 传 热 


160 
] —s— Normal 
1404 —e— MSLB 
S ] —4— FLOW 
E 1204 
$ 1004 
$ 80 
z 
3 60 
> ] 
404 
0 5 10 15 20 25 30 


Flow direction (axial direction) / mm 


11 又 排 模型 外 壁面 轴 向 应 力 分 布 特性 


Fig.11 Characteristics of axial stress distribution on 


the outer wall of cross model 


进行 了 正常 运行 和 典型 事故 参数 下 的 结构 响应 计算 和 分 析 ， 形 成 的 主要 结论 如 下 : 


1) 


壁面 ， 在 一 、 二 次 合流 动 条 件 相同 的 情况 下 ， 又 排 模型 在 正常 


AUREUS 


排 模型 ; 


2) 


3) 


壁面 的 平台 应 力 值 


管 排 列 方式 模型 ， 


等 效应 力 在 轴 向 近 开 似 呈 以 中 心 模 截面 的 对 称 分 布 ， 内 壁面 应 力 水 平整 体高 于 外 
运行 和 典型 事故 参数 下 的 应 力 和 变形 大 于 顺 


传 热管 根部 横 截面 的 等 效应 力 沿 壁 厚 方向 先 减 小 再 增 大 ,应 力 最 大 值 在 内 壁面 ， 最 小 值 在 壁 厚 中 部 
靠近 外 壁面 位 置 ， 传 热管 中 部 横 截面 的 等 效应 力 沿 壁 厚 方向 呈现 近似 线性 降低 ， 内 壁面 应 力 最 大 ; 
传 热管 根部 区 域内 、 外 壁面 的 等 效应 力 沿 轴 向 高 度 先 迅 速 降低 而 后 快速 升 高 ， 并 形成 应 力 平台 ; 内 
低 于 根部 应 力 ， 而 外 壁面 的 平台 应 力 值 高 于 根部 应 力 。 
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文章 修改 ; 朱 光 昱 负责 调 下 


及 文章 校对 ; 唐 晖 晖 负责 


参考 文献 


1 ” 王 莹 杰 , 王 明 军 , 王 树 强 ,等 


技术 指导 ; 刘 巧 凤 人 负责 


文章 校对 。 


. 核 动力 系统 蒸汽 发 生 器 U 形 管 流 致 振动 数值 模拟 [四 . 核 科学 与 


[ 程 ,2022,42(03):723-730. 


WANG Yingjie, WANG Mingjun, WANG Shuqiang, et al. Numerical Simulation of Flow-induced Vibration of U-tube in Steam 
Generator of Nuclear Power System [J]. Nuclear Science and Engineering,2022,42(03):723-730. 


2 WEE, P 


3 MEF.ERAEŠ 


清华 


pK. 1% 


电厂 系统 及 设备 [M]. 北京: 


> 


华 大 学 出 版 社 ， 
ZANG Xinian, SHEN Shifei. Nuclear power plant systems and equipment[M]. Beijing:Tsinghua University Press, 2003:65-84 


TCRUESIHZRTAUR ERFASST 43 CJ] 422) 7] 1:86,2000,(02):152-156. 


2003: 65-84 


ZANG Xinian. Event Tree Analyses of Steam Generator Heating Surface Tube Rupture Induced by Main Steam LineBreak[J]. 
Nuclear Power Engineering,2000,(02):152-156. 


4 ”张骞 . AP1000 2&1 AE TE HAVE (SGTR) SE LAT 
ZHANG Qian. 7| 


Sha 


nghai: Shanghai Jiao Tong University,2015. 


5 BH 


并 轩 , 张 竞 宇 ,王晓东 ,等 . 压 水 堆 核 


AD]. 上 海 : 上 海 交 通 
KÆ. STUDY ON THE STEAM GENERATOR TUBE RUPTURE ACCIDENT OF AP1000 NPP AP1000[D]. 


学 ,2015. 


电厂 蒸汽 发 生 器 传 热管 道破 裂 事故 源 项 的 计算 分 析 四 . 核 技 术 ,2020,43(06):31 


-36. 


FAN Yuxuan, ZHANG Jingyu, WANG Xiaodong, et al. Calculation and analysis of steam generation tube rupture accident 
source term in PWR [J]. Nuclear Techniques,2020,43(06):31-36. 


XXXXXX-7 


HE dX 术 20XX,XX: XXXXXX 


10 


11 


12 


KUEHLERT K, WEBB S, SCHOWALTER D, etal. Simulation of the fluid-structure-interaction of steam generator tubes and 
bluff bodies[J]. Nuclear Engineering and Design, 2008, 238(8): 2048—2054. 

Th 2 ARE RA Sb SCR SE. AE F Uic ARA OI ZRT E A RES EJI 4 BP LT] £5 T 5538,2014,65(S1):364-370. 

SUN Baozhi, ZHENG Lusong, HAN Wenjing, et al. Analysis on structural stress of tube in steam generator based on 
]. Journal of Chemical Industry and Engineering,2014,65(81):364-370. 

郑 陆 松 , 孙 宝 芝 , 杨 元 龙 ,等 .基于 流 热 固 耦 合 的 核电 蒸汽 发 生 器 传 热管 热 应 力 数 值 模拟 四 .原子 能 科学 技 
术 ,2014,48(01):74-80. 

ZHENG Lusong, SUN Baozhi, YANG Yuanlong, et al. Numerical Simulation on Thermal Stress of Tube in Steam Generator of 


fluid-structure interaction[J 


È 


Nuclear Power Plant Based on Fluid - thermal - structure Interaction [J]. Atomic Energy Science and 
Technology,2014,48(01):74-80. 

RER, He, ] Ade, SE 3E T UA DSL GRO U IHESU S P P 5 IL 739 2 TET] EGCS,2023,46(01):109-116. 
CHEN Yushuang, TIAN Jian, DING Mengting, et al. Temperature field and stress field analysis of molten salt U-tube heat 
exchanger based on fluid-thermal-solid coupling method [J]. 213 75,2023,46(01):109-116. 


Sun B Zheng L , Yang L , et al.A co upled stress analysis of the steam generator tube considering the influence of the fluid flow 
and heat transfer in the primary and secondary sides[J].Applied Thermal Engineering,2015,87803-815. 

BIAR, IEZ, SE UIT ER TAE S P AVE A I V [8 e UE DA M 548,2016,67(81):217-223. 

ZHAO Yingjie, SUN Baozhi, SHI Jianxin, et al. Two-way fluid solid interaction numerical analysis of stearm generator heat 
transfer tube[J]. Journal of Chemical Industry and Engineering,2016,67(S1):217-223. 

ALIA E STAR AE ti Pe AVES RO I2 7J 4T [D ]. 636 r8 7] $2015. 

JIANG Shuang. Stress Analysis of Heat Transfer Tube for PWR Steam Generator [D]. North China Electric Power 
niversity,2015. 


IIT 


[ani 


